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1 Introduction 

In 1962, Josephson £Q predicted a novel phenomenon in superconductivity, 
namely when two different superconductors were brought into close contact. 
Based on elementary quantum mechanics, he predicted the existence of a su- 
percurrent with a peculiar current-voltage dependence. He argued that there 
would emerge a current of Cooper pairs which is proportional to the sine of 
the phase difference of the order parameters of both superconductors. The 
success of this prediction was immediate when indeed this phenomenon was 
experimentally observed one year later It counts as one of the greatest 
successes of quantum mechanics in physics and you will find a chapter on 
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the Josephson effects in almost every textbook on superconductivity. The 
increase of knowledge on this subject in theoretical solid state physics in the 
following decades has been tremendous and applications of Josephson junc- 
tions in electronic devices have been developed [3]. Progress in conceiving a 
microscopic theory for the Josephson effect in rigorous quantum statistical 
mechanics was made when Sewell obtained the Josephson and Meissner ef- 
fects in an model independent approach from the assumption of off-diagonal 
long range order and local gauge covariance [I] . 

In jSj we considered a concrete microscopic quantum model yielding an 
ab initio and rigorous understanding of the emerging of a Josephson current, 
which is computed and numerically calculated. The model consists of two 
two-dimensional superconducting plates having a common one-dimensional 
contact surface through which Cooper pairs can tunnel and as such induce a 
current. In section 2 we repeat the essentials of the model for which we con- 
struct a non-equilibrium steady state (NESS). One of the attractive aspects 
of the construction is that our NESS has the nice property of having a finite 
interaction area. We derive an analytical expression for the current in the 
case that the phases of the two superconductors are not too large. We find 
back the perfect sine-behavior. 

Section 3 is devoted to the study of the appearance of Goldstone bosons 
in the junction, due to the interaction of the two bulk superconductors and 
a direct consequence of the gauge symmetry breaking. We apply the general 
result of [Zj where one finds the general explicit construction of the normal 
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Figure 1: Two superconductors wit a contact surface 



coordinates of the Goldstone boson as a consequence of the spontaneous 
symmetry breaking which appears in our model. The Goldstone bosons of 
the bulk superconductors can also be found in that paper. Here we find two 
supplementary Goldstone bosons. We construct their normal coordinates and 
their dynamics induced by the micro-dynamics of the model. We consider 
their dynamics in diagonal form and analyze its spectrum, again as in section 
2 as a function of the phase difference of the two bulk superconductors. Here 
one finds a cosine-behavior. 



2 Micro model for Josephson Junctions 

As said above we consider the model jH] for two two-dimensional supercon- 
ducting plates / and II with a common one-dimensional contact surface (line) 
through which the Cooper pairs can travel in order to induce a current. 



The two superconductors are modeled by the strong coupling BCS-model 
on a square lattice using the Anderson quasi-spin formalism and described 
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by the Hamiltonians Hi ^ with % — I, II 

N 1 N 

H L N = ^o, o - Yl l ) a i( m > n )' e * > 

k,l=l k,l,m,n=l 

acting on the Hilbert space (Ei^C^, of and a z = crfcrj — af are copies of 
the Pauli matrices. The a + and a~ represent the creation and annihilation 
operators of the Cooper pairs of the superconductors; the e do represent the 
kinetic energies of the Cooper pairs. 

The junction between the superconductors I and II is modeled by the 
interaction 

N 

Vn = -JJ 1)^7/(^2,1) + h.c), 7 >0 (2.2) 

fci,fe 2 =i 

which is responsable for the Cooper pair tunneling through the barrier. A 
pair at the site (kj, 1) of the first superconductor can tunnel through the 
junction and create a pair at the site (k 2 , 1) of the second superconductor 
and vice- versa. The coupling constant 7 governs the rate of this process. Note 
that only Cooper pairs on the contact surfaces of / and II participate in this 
process. Remark that only N sites of each superconductor are interacting 
with N sites of the other one. The lattice permutation invariance of the 
Hamiltonians ()2.1|) and the interaction ()2.2j) make the model given by the 
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total Hamiltonian of the system 



H N = H IN + H II>N + V N (2.3) 
exactly soluble in the thermodynamic limit N tending to infinity 

2.1 Equilibrium states of the non- interacting super- 
conductors 

For completeness we discuss here the equilibrium states of the Hamiltonians 
(|2.1jl . We treat the first (I) one explicitly, the second is analogous and ob- 
tained by replacing the index / by the index II. Following [B], the extremal 
equilibrium states at inverse temperature fli in the thermodynamic limit are 
the product states ui tpi with the expectation values of all tensor product ob- 
servables X = X X1 ® X X2 <g) Xi,x 2 , ... G N 2 and all X Xj G M 2 (2 by 2 
complex matrices), given by 

m^o= n Tr ^-^ ( 2 - 4 ) 

Here p ipix is the x-copy of the 2 by 2 density matrix p Lpi G M 2 , solution of 
the selfconsistency equation 



with the one-site effective Hamiltonian h ipi given by 

h ipi = e I (T z I -\ I (e i * I (Tj + h.c.), A 7 > (2.6) 

Clearly the density matrix l2*31 is also determined by the equivalent selfcon- 
sistency equation for the order parameter Xj = \u ipi (ay)\ which by explicit 
computation becomes 

A/(l - — tanh pjkj) = 0, m= ^ej + Xj (2.7) 

Remark that {±/x/} constitutes the spectrum of the effective Hamiltonian 
h ipi which is independent of the phase angle ipj. 

It can readily be seen that 12. II admits always a solution Xj = 0. It yields the 
normal phase state of the superconductor. For £j < | and (3i large enough 
or temperature low enough, there exists also a solution A/ ^ 0. These so- 
lutions yield the superconducting phase states. The phase ipj can be fixed 
freely, yielding an infinite degeneration of the equilibrium states under the 
conditions mentioned above. 

The second superconductor(II) has analogous phase states with phases which 
can be chosen independently from the first one. 

In the following we fix for each superconductor such a superconducting phase 
state denoted by Wj m and 0Ju^ n with <pj ^ ip n . 
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2.2 Non-equilibrium steady state (NESS) 

Now we construct a non-equilibrium but steady state (NESS) for the total 
interacting system 12.31 We start from the product state uj = uj[ jipi ® wn,<p u 
on the system, i.e. the state of the system of the two superconductors in 
their respective superconducting phase states characterized by their phases, 
inverse temperatures and kinetic energies. In this state uj one can compute 
the global dynamics yielding a time evolution 

a t (.) = uj — lim(exp UH N . exp itH N ) (2.8) 

where uj — lim is the weak limit under the state uj. Now we are ready to look 
for a state uj which is invariant under the dynamics 13.131 
Due to the specific lattice permutation symmetry of the model, it is natural 
to choose this state uj among the product states i.e. 

Q VI (X)=Tr J] p x X x (2.9) 
zeN 2 

As the state should be time invariant, it should have the same lattice invari- 
ance as the Hamiltonian 12.31 The symmetry of the Hamiltonian divides the 
total system into four parts. 

There are the bulk parts of the superconductors I and II which we denote 
by I a and II a , and there are the contact or surface parts denoted by lb and 
I If,. Therefore one can write the state uj as a tensor product of four symmetric 
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Figure 2: Division of the system in four subsystems 

product states on their different regions: 

oo = u Ia ® tu Ib ® u IIb ® a> 7a (2. 10) 

Finally we require that the state w is a steady state 

\imu)([H N ,X\) = (2.11) 

Due to the product structure of the state, the Hamiltonian can be identified 
with an effective Hamiltonian [Oj of the type Hn = J2 x ei N uii N ^ x > where 
h x G M 2tX . Imposing this time invariance yields 



hi = < 



e 7 a z (z) - u{aJJa Ia (t) - uj{a la )aj a {t), 



(2.12) 



tu° z u h {%) ~ u{aJ Ia )a IIb {i) - co(a H JaJ Ib (i) 
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We use also the notation 



A/ = Xiexpiip! = u(aj b ), y?/ = argu)(cr+) 



(2.13) 



and analogously for the second superconductor with / replaced by II. 
The local density matrices p x for x G / U II of the state u) are the h x invariant 
projections of the density matrices p x of u (12 .4j) . For more details of the 
construction, see [5]. In any case, for x G I a U // a : p x = Px as follows 
from 12.121 For the lattice point x G /> U //ft one readily computes the 
selfconsistency non-equilibrium equations 



Together with the selfconsistency equations 12. ll for Xj and Xjj, the equations 
12 .141 form a set of six coupled equations whose solutions determine the non- 
equilibrium steady state u of the total system. 

This state divides the system into four parts. The bulk parts of both super- 
conductors away from the contact surface do not feel each other nor do they 
feel the surface. They behave as stable reservoirs. On the contact surfaces Ij, 
and lib the system is effectively perturbed and influenced by the properties 
of the states of both superconductors. 

Remark that we limited the interaction to take place only on a contact sur- 
face of one layer thickness. It is clear that the whole construction can be 




ej + AjjA/e^ + jX n e i,f>II \ cos(0j - fa) 

ej + |A/e^ + 'jX II e i ^'\ 2 
) eh + XnlXne^' + 7 A 7 e^ | cos(0 J7 - 4> H ) 
e 2 n + |A//e^« + 7 A/e*| 2 

(2.14) 
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generalized to the case of any fixed finite number of layers. 
In jHj we considered the currents of Cooper pairs emerging in the system 
traveling from the superconductor I to II and vice versa. One considers the 
relative particle number operator of the Cooper pairs 

Qn= J2^7-^7i) (2-15) 

X6JV 2 

The local relative current is then 

2' N 

J(Q N ) = i[H N , Q N ] = - J2 5>7 6 (i, l)a+ b {j, 1) - h.c.) (2.16) 

Remark that in the observable current there is no direct contribution 
from the bulk of the two superconductors, only the two contact layers are 
contributing. Moreover one remarks that this current is of the same order of 
magnitude as the contact surface, namely N. 

The Josephson current measured in the thermodynamic limit (N — > oo) state 
oj, called NESS, is readily calculated and given by 

j(Q) =lim Uj{J ^ N)) = -4 7 A/A// sin - (p n ) (2.17) 

This result was obtained in [Sj, here it is written in a more concise form. 
However one has to realize that the quantities A and (p are functions of the 
originally given parameters X I ,X II ,(p I and ipu, given bv !2.14l First of all it is 
easy to see that these equations are shift invariant for an arbitrary shift of the 
two originally given angles. This means that, without loss of generality, one 
can take one of the angels, say fi, equal to zero. Furthermore, it is natural 
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to assume the coupling constant 7 to be very small, i.e. 7 <C min(ej, en). 
Therefore it is reasonable to compute the quantities A and (p only up to first 
order in this parameter 7. 

By multiplying the two selfconsistency equations with each other, using the 
fact that tpj = and taking tpjj < | one gets that the phase difference 
Pi — 'fill is proportional to the phase difference tpj — tp n . Suppose now that 
also the second phase tpjj is small, corresponding to the usual experimental 
regime. Then one remarks that also tpj = 0. It follows that 



ip~ H « tp u (2.18) 
After substitution of 12. 181 in 12.141 one gets 

A / = A / -7^ (2.19) 

\ \ I A I e II 
Mi = A// + 7 — T~ 

V11 

After substitution of all these equations in the formula 12.171 one gets the 
expected formula for the Josephson current 



j{Q) = -^XAn sh% 7 - cp n ) (2.20) 

yielding an analytical expression for the current for small phase differences 
between the two bulk superconductors. A numerical computation of the 
current for arbitrary phase differences is found in 
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3 Symmetry breaking and Goldstone bosons 

As is well known, spontaneous symmetry breakdown (SSB) is one of the 
basic features accompanying collective phenomena. It became a representa- 
tive tool for the analysis of many phenomena in modern physics. For long 
range interactions, it is typical that SSB is accompanied also by the breaking 
of the symmetry of the dynamics. The latter phenomenon is known to be 
accompanied by the occurrence of oscillations of a Goldstone boson with a 
non-vanishing energy spectrum. These oscillations together with the Gold- 
stone boson disappear if the SSB disappears. In [7j one was able to construct 
explicitly the normal coordinates of these new particles called Goldstone par- 
ticles. In particular for mean field systems such as the BCS-model jB], the 
Overhauser model [5], a spin density wave model [TO], the anharmonic crystal 
model [in], and the jellium model [TT| . one has constructed these Goldstone 
boson normal coordinates. 

Our two-dimensional model consisting of two interacting superconductors 
also shows the phenomenon of SSB. As the main contribution of this paper 
we consider the construction and the calculation of the spectrum of the cor- 
responding Goldstone bosons. 

As far as the Josephson current, computed in the previous section, is con- 
cerned we remark that the bulk parts I a and II a of the superconductors are 
not contributing to it. Therefore it is reasonable to look for the Goldstone 
particles within the contact areas and lib- In particular we compute the 
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normal coordinates and the dynamics of the Goldstone bosons of this junc- 
tion area. From [2j we know that the canonical coordinates of these bosons 
are given by the fluctuation operators of the generator of the broken symme- 
try and of the order parameter operator. 

As the following gauge transformation holds 

e iaaZ a + e- iaa * = a + e 2ia ; a G R (3.1) 

the a z are the local generators of the broken gauge symmetry of the effective 
Hamiltonian 12. 12l Indeed for all a 

Cj(e iaaZ a + e- iaa ~~) = Co{a + )e 2ia (3.2) 

proving that the state is not invariant under the U{1) gauge group. 
Therefore we consider the local operators, for i G lb, Ih 



Qi = ^4fW + S> (A^ + + h.c.) (3.3) 



where for all i G h' 



Ai = A Ib = u(aj) + iu{a IIh ) (3.4) 



/2 i = // Ii = v 'e I 2 + |A J J 2 (3.5) 
e/ = Cj (3.6) 

and equivalently with i G lib i-e. by substitution of / by II and vice versa. 
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Remark that the operator Qi is indeed the generator of the gauge trans- 
formations, namely up to a constant equal to a z , but normalized to zero 
expectation value 



We consider also essentially the order parameter operator a fluctuation 

P s = f (A>+ - h.c.) (3.8) 

Again remark that d>{Pj) = i.e. also this operator is duly normalized to 
zero. 

Using the general quantum fluctuation theory for product states [12]; one 
computes the following quantum central limits in the given state u and obtain 
the normal coordinates of two Goldstone bosons. For the region lb one gets 
the normal coordinates 




^Co(a z ) + 4^2|A| 2 = 



(3.7) 




(3.9) 




and for the region Ii& one gets the normal coordinates 



bn b (Q) = lim 



(3.10) 



b IIb (P) =lim^^E. 
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In 13 .91 one gets the normal coordinates of a first Goldstone boson, and in 13. 101 
the normal coordinates of a second independent Goldstone boson. Indeed, 
by a straightforward computation one checks readily the following canonical 
commutation relations 

\b Ib {Q)Mi b {Q)] = [hMM h {P)]= (3.H) 

[b Ib (P),b Ih (Q)} = [bi b (P),bn b (P)] = 

[b n SQ)M h (P)] = ^jt b 

Remark that in the case of temperatures above the critical ones of the bulk 
superconductors the order parameters vanish: Xj b = Xn b = such that all 
commutators in 13.111 vanish. Also one has 

u{b Ib {Qf) = u{b IIb {Q) 2 ) = u>{b h {P?) = u{b Ilb {P?) = (3.12) 

and hence all the operators themselves vanish: bj b (Q) = bj b (P) = bn b (Q) = 
bn b {Q) = 0, i.e. the Goldstone bosons disappear in the normal phases. 
Next we consider the dynamics of the Goldstone bosons in the case of su- 
perconducting phases for the bulk superconductors. We consider the time 
evolution of the normal modes 13.31 and 13.81 which is induced by the initial 
micro-dynamics given by the effective Hamiltonian 12. 121 In general, let A be 
a local observable situated at the lattice point x G h or II b , then denote d t 
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the time evolution of the fluctuation of A after time t. It is given by 

a t b h (A) = b Ib (exp(ith x )Aex.p(-ith x )) (3.13) 
a t b IIb (A) = b IIb (exp(ith x )Aexp(-ith x )) 

Of course the operator A stands for the operators Q x ()3.3|) and P x ()3.8|) . 
A straightforward computation of the dynamics using ITTTTTTl yields the simple 
solutions 

a t b Ib (Q) = b h (Q) cos(2/2 4 i) + b h (P) sin(2/i /6 t) (3. 14) 

a t b h (P) = -b Ib (Q) sm(2fi Ib t) + b h {P) cos(2/2/ 6 t) 

and analogously for the surface one gets the same dynamics for the second 
Goldstone boson by replacing the index If, by the index 11^. 
The two bosons behave dynamically as two independent quantum harmonic 
oscillators with frequencies 

f>i = 2/l Ib (3.15) 
i>u = 2p>ii b 

For the bulk superconductors I a and II a the Goldstone bosons dynamics 
were computed before jH]. The frequencies Uj = 2/ij and v u = 2fi n are 
clearly phase independent. 

However for the frequencies of the Goldstone bosons considered in this paper, 
the situation is completely different. The frequencies 13.151 computed above 
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do depend on the phase difference tpi — (fn of the two phases of the bulk 
superconductors. 

We consider the phase dependence explicitly for the region 7?,, the com- 
putation for the second region is analogous: 



Remark that the A/ and (pn are determined by the parameters ej, e/j, A/, 
Xij and the phases cpi and <pn through the selfconsistency equations 12.141 

It is instructive again to get an explicit form of the frequencies in terms 
of the given parameters of the system. We derive again the formula for the 
frequency vi up to first order in the coupling constant 7, and in the case 
that the phase difference tpi — cpu is small, a situation explored before for 
the current. 



Hence one gets the following expressions for the dynamical frequencies of the 
Goldstone modes 




(3.16) 



We get from l2~T£l l2~T31 and l3~T51 



fx 2 ! = ej + A^ + 27A/A// cos(y?/ - (p n ) 



(3.17) 



Vl = Vi + 47 




cos(^/ - ipii) 



cos(v? 7 - ip n ) 



(3.18) 



(3.19) 
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From these expressions the dependence of the Goldstone frequencies on the 
phase differences of the bulk superconductors is explicitly given. We learn 
that these frequencies decrease when the phase difference increase in con- 
tradistinction with the current. The current has a sine-behavior, the frequen- 
cies a cosine-behavior. This point may be interesting from the experimental 
point of view. 
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